Transport of a heated granular gas in a washboard potential 
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We study numerically the motion of a one dimensional array of Brownian particles in a washboard 
potential, driven by an external stochastic force and interacting via short range repulsive forces. In 
particular, we investigate the role of instantaneous elastic and inelastic collisions on the system 
dynamics and transport. The system displays a locked regime, where particles may move only via 
activated processes and a running regime where particles drift along the direction of the applied 
field. By tuning the value of the friction parameter controlling the Brownian motion we explore 
both the overdamped dynamics and the underdamped dynamics. In the two regimes we considered 
the mobility and the diffusivity of the system as functions of the tilt and other relevant control 
parameters such as, coefficient of restitution, particle size and total number of particles. We find 
that, while in the overdamped regime, the results for the interacting systems present similarities 
with the known non-interacting case, in the underdamped regime, the inelastic collisions determine 
a rich variety of behaviors among which is an unexpected enhancement of the inelastic diffusion. 

PACS numbers: 05.40.-a,61.20.Gy.05.10.Gg 



I. INTRODUCTION 

A large class of phenomena in biology, chemistry, en- 
gineering and physics occurs via the transport of par- 
ticles driven along periodic substrates by an external 
bias. These phenomena include polymer diffusion at 
interfaces^, motion of fluxons in superconductors con- 
taining a periodic arrangement of defects -^^4^ adsorption 
on crystal surfaces^, super-ionic conduction^, motion of 
molecular motors along microtubules^, granular flows on 
rough inclined substrates— i^i^i^iii or in a narrow pipe^^. 

The study of transport properties of granular systems 
represents an open issue in Statistical Mechanics of con- 
siderable interest and difficulty in view of the continuous 
energy dissipation caused by particle inelastic interac- 
tionSii^. In this perspective, we consider a simple model 
consisting of a granular fluid moving on a tilted rough 
substrate, which may favor clustering and jamming be- 
haviors. 

In our formulation, a granular system is a large number 
of particles (grains) colliding with one another and losing 
a little energy in each coUisio n^^i^^i^^i^^i^^ . If such a sys- 
tem is shaken to keep it in motion its dynamics resembles 
that of fluids as the grains move randomlyii^ 

We carry out a comparative study of the behavior of 
three models: the inelastic particle system (IPS), the 
elastic particle system (EPS), and flnally the non inter- 
acting system or independent Brownian particles (IBP) 
to understand how the interactions influence the collec- 
tive transport. The absence of interactions between par- 
ticles makes the basic phenomenology of the IBP well 
understood. It reduces, indeed, to the motion of a sin- 
gle particle in a force field resulting from the interplay 
between a constant force F , such as gravity or electric 
field, a spatially periodic force, simulating the presence 
of a rough substrate, a viscous force accounting for the 



friction, and a noise term, representing thermal fluctua- 
tions The scenario is the following: particles diffuse 
with a bias in the direction of the steady force, F, with 
an average velocity, which is an increasing function 
of the applied load F. However, when the load is be- 
low a critical value, Fc, and the noise sufficiently small, 
the average velocity tends to zero. In other words the 
particles remain locked in the minima of the periodic po- 
tential. On the other hand, above the threshold, F^ 
the particles may travel from one minimum to the other 
along the tilt direction. The critical value, Fc, depends 
on roughness, friction and temperature. According to 
the value of the friction coefficient two different scenar- 
ios may be observed. In the so called overdamped limit, 
Fc occurs only when the tilt is such that the local minima 
completely disappear. In the opposite limit, i.e., in the 
underdamped regime, the particles may overcome many 
barriers even for low values of F at which the potential 
still displays local minima. In fact, the particles may 
cross a barrier, provided their gain in potential energy, 
to go from one minimum to the next, exceeds the energy 
dissipated along this pathway. In the underdamped limit 
the exact value of Fc depends on the temperature. 

Two indicators may serve to characterize the collective 
transport properties of the system, the mobility, /i, and 
the diffusion constant, D. The former relates the average 
velocity to the tilt F according to (w) = fi{F)F. The 
latter measures the average spreading of particles: 



R{t) ^ {[xit) ~ (xcAiit))? 



(1) 



The average is meant over the thermal noise realizations. 
In all the cases where R{t) asymptotically grows linearly 
in time, we can identify the normal diffusion constant D 
defined as the slope of the law 

R{t) - 2Dt 
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that generally is computed numerically through simula- 
tions of the system evolution. The study of the depen- 
dence of fi and D on F and other system parameter is 
necessary to determine the efficiency of the transport, 
i.e. the ratio between drift and spreading of a swarm of 
particles. 

A detailed theory capable of describing the behavior 
of the IBP for arbitrary values of the damping has been 
developed by Risken2£i2i who derived an analytical ex- 
pression for the mobility, /i, as a function of noise level, 
tilting force and substrate characteristics. A closed for- 
mula for the diffusion coefficient exists only for the over- 
damped situation^^. It predicts the presence of a striking 
enhancement of D near a threshold, F3, separating the 
locked from the running phase and corresponding to the 
disappearance of the local minima of the potential. 

On the other hand, in the underdamped regime there is 
no analytic expression for D, but the enhancement of the 
diffusion coefficient, occurs at a lower threshold, F^^'^^ 
due to the effect of inertia, and the phenomenon remains 
qualitatively similar. 

The goal of the present paper is to consider the effects 
of interactions among particles on the transport proper- 
ties in such systems. Specifically, we will address the fol- 
lowing question: how do interactions occurring via elas- 
tic/inelastic hard-core collisions affect the mobility and 
the diffusion of an assembly of particles? This issue is 
particular relevant in one dimension where a repulsive 
hard-core interaction inhibits the particles to pass each 
other; this constraint is known to inffuence dramatically 
the dynamics of a group of particles. As an example, 
we can mention the anomalous self-diffusion in single file 
systems-^ 

We focus our dynamical approach only on exploring 
the effects of repulsive interactions via impulsive contact 
collisions. We consider two possibilities, energy conserv- 
ing collisions and dissipative inelastic collisions. Both sit- 
uations are still largely unexplored in systems with wash- 
board potentials and display non-trivial behaviors, as we 
shall illustrate. For instance, the mutual repulsion be- 
tween the particles induces dynamical correlations which 
may favor or hinder their motions: it can promote their 
exit from a potential well via energetic collisions, or, on 
the contrary, it can prevent a jump to a given site when 
this is too crowded. In addition, the granular tempera- 
ture, defined as the average kinetic energy per particle, in 
inelastic systems is in general lower than the correspond- 
ing temperature of elastic systems. Therefore, the trans- 
port, occurring via thermally activated processes across 
substrate barriers, is expected to be less efficient for in- 
elastic particles^l. Then few questions can be addressed. 
Does the diffusivity D present an enhancement analogous 
to that of non interacting systems? Does the threshold 
locked-to-running occur at larger values of the tilt? In 
the following, we answer the questions by examining a 
variety of situations and analyzing the overdamped and 
the underdamped regimes separately. 

The paper is organized as follows. In section II, 



we present the model and discuss the contact elas- 
tic/inelastic interaction with their main implication in 
the system dynamics. In Section III, we analyze the re- 
sults from the molecular dynamics (MD) simulations of 
the over-damped regime. In section IV, we summarize 
the corresponding results for the under-damped regime 
discussing the salient differences with the over-damped 
case. Finally, conclusions are drawn in section V. 



II. MODEL 

We consider a randomly driven granular gas, already 
introduced in previous works^^, composed of N impen- 
etrable hard-rods, of mass m and size a, moving on a 
line of length L in the presence of an external washboard 
potential. The overall dynamics of this gas is described 
by the Langevin equation for each rod 



d?Xi{t) dxi{t) 



= —7717- 



dt 



dxj 



that embodies four types of physical phenomena: friction 
with the surroundings, random accelerations due to ex- 
ternal driving, inelastic collisions among the particles and 
external time independent, but spatially varying, force. 
We model the first two effects by means of a stochastic 
bath with a viscous friction —m'jii and Gaussian random 
force, with zero average and covariance 



(3) 



satisfying a fluctuation-dissipation relation, with T pro- 
portional to the intensity of the stochastic drivingSS. The 
damping term fulfilling Einstein's relation renders the 
system stationary even in the absence of collisional dissi- 
pation and physically can represent the friction between 
the particles and the container. The interactions among 
the particles, indicated formally by J^j fij iii Eq. ([2]), 
amount to simple impulsive forces acting on the rod i due 
to the collisions with the neighboring particles j30,3i,32,33 ^ 
Thus a rod i performs, under the influence of the bath 
and collisions only, independent Brownian trajectories, 
unless it gets in contact with particle j, i.e. \xj —Xi\ = cr, 
at which point the velocities of the colliding pair 
change instantaneously according to the inelastic rule 
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the prime indicating post-coUisional variables and a the 
coefficient of restitution, with < a < 1. 

The total force, in Eq. on each particle is supple- 
mented by terms representing the presence of a constant 
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external bias, F, and a periodic substrate generating a 
potential ^'(x), with period w and amplitude Vq, 



'^{x) = — Vb cos(27ra;/w). 



(4) 



A quantitative criterion to distinguish between an over- 
damped and an underdamped regime is provided by the 
dimensionless parameter F — 7/cjo, where 



uJo = \j — 
m 



dx"^ 



is the oscillation frequency at the potential minimum. If 
r ^ 1 the dynamics is overdamped, and in the opposite 
limit is underdamped. 

As discussed by Borromeo and Marchesoni^^^, the 
present model is equivalent to consider the particles in 
a traveling potential, i.e. a periodic potential moving 
at speed c, which is related to the tilt force by the re- 
lation F — mjc. Interestingly for speeds lower than a 
certain threshold the traveling wave has the capability 
of dragging the particles, a mechanism known as Stokes' 
drift^*^''^^ . 

Thus, the system can equivalently describe the physics 
of a one dimensional granular gas where the periodic po- 
tential 'I'(x) represents a series of compartments sepa- 
rated by walls, an experimental set-up recently employed 
to study the clustering behavior of vibrated granular 
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or the roughness of an inclined planeSi 



8^ 



III. NEARLY OVERDAMPED REGIME 

We shall begin by considering the nearly overdamped 
regime (F ~ 2.1), whose study is better understood be- 
cause the system reaches a steady state rapidly due to 
the large value of the friction. 

Our MD simulations were carried out by evolving an 
initial configuration, where the particles were all located 
in the central well without overlaps and their velocities 
were extracted from a Maxwell-Boltzmann distribution 
of temperature T. Each run involved N = 256 particles 
of size a = 1, mass m = 1, with three different values 
of the coefficient of restitution, a = 0.8, 0.9, and 1.0, 
respectively. The substrate potential was characterized 
by wells of width w = 40Qa and height Vq = 9.0, in units 
of fcsT. Finally, the heat-bath temperature has been 
chosen to be T = 1. and fc^ = 1. The time is measured 
in units of i„ = cry^m/T. The friction coefficient is 7 = 
2/i„. The cyclic system studied has been taken to be of 
length L = 10^ w so that it is virtually equivalent to a 
system with open boundary conditions. In this condition 
of high dilution, the global system density p — N/L is 
extremely low, however it is not meaningful parameter, 
rather it is the initial density profile, characterized by 
the number of particles in a well pw — N/w, that has 
a strong influence on the system behaviour. Indeed, at 
the beginning of the evolution, the system needs to be 
packed enough to reach a not negligible collision rate. 



Only in the later stages of the simulations, a crossover is 
observed toward the behaviour of rarefied gases. 

At first sight, the MD simulations of the interacting 
system display a mobility quite similar to that of the 
non-interacting system: the curves 7/i versus F for elas- 
tic and inelastic hard-rod systems are shown in Fig. [TJ In 
the same figure, we also plot the corresponding quantity 
relative to the IBP as calculated numerically by means 
of the continued fraction method'^°. However, a closer 




FIG. 1: Dimensionless mobility 7/i as a function of the 
rescaled external force F/F3 {F3 = 2nVo/w) for a — 1.0 
(circles) and a = 0.8 (squares) for a system with A*' = 256 
particles, temperature T = 1, friction 7 = 2/tu, and potential 
amplitude Vb = 9.0, cr — 1, w/a = 400. The system evolution 
is simulated for t,nax ~ 10000 time units. The full line indi- 
cates the corresponding curve for the non interacting particles 
obtained by the continued fraction method^. 

inspection reveals that differences do exist between the 
IBP and the interacting systems with different inelas- 
ticity a. Figure [2] shows the ratio of the mobilities of 
interacting systems to the mobility of the IBP. This ra- 
tio, for small values of the load, F <^ F3 — 27rVb/w, can 
be significantly different from 1, and more specifically 
the mobility of the EPS exceeds the mobility of the IBP, 
whereas the mobility of the IPS is lower. The reason for 
these differences can be found in the fact that contact 
interactions may change drastically, with respect to the 
IBP, the time that particles spend in a given potential 
minimum. 

In more detail, these behaviors can be explained by 
recalling that two competing mechanisms contribute to 
the system mobility: the excluded volume and the in- 
elasticity. The first leads to an effective reduction of the 
barrier height and favors the escape from the wells, thus 
basically increasing the mobility. The second, instead, 
tends to decrease the average kinetic energy, rendering 
longer, on average, the time spent by the particles in a 
given well. For low values of the forcing field F, the mu- 
tual repulsion dominates and thus we observe a larger 
mobility of the EPS with respect to the IBP at the same 
values of F and T. However, as we switch the inelasticity 
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FIG. 2: The ratio between mobilities of interacting and non 
interacting systems as a function of the rescaled external force 
F/F3 {F3 = 2-n-Vo/w) for a = 1.0 (circles), a = 0.9 (squares) 
and a — 0.8 (triangles). The remaining parameters are the 
same as in Fig. [1] 



on, the mobility of the IPS can become lower than that 
corresponding to the IBP. For F > F3, /i turns out to be 
the same for all systems because in this regime only the 
forcing field F matters. In fact, the absence of minima 
in the total potential reduces the influence of the mutual 
interactions so the particles follow coherently the strong 
effect of the drift. 

The role of the excluded volume can be highlighted 
by monitoring the trajectories of some specific tagged 
particles whose dynamics can be very representative of 
the evolution of the whole system. More specifically, we 
consider the two extremal particles, 1 and N (having 
labeled the rods along the drift direction from 1 to N) 
and the central particle i = N/2. 

The first particle is, on average, braked by the colli- 
sions of the preceding particles, whereas the A^-th particle 
is pushed ahead by the pressure of those behind. Their 
overall behavior is expected to be qualitatively rather 
different from that of the central particle. In figure [3l 
we display, for a = 1.0, 0.9 and 0.8, the evolution of the 
mean displacement of these three tracer particles from 
their initial position. Each displacement is normalized 
with respect to the displacement of the center of mass of 
the whole system from its initial position: 



Zk{t) 



(|xfc(t)-xfc(0)|) 

{\xcM{t)-XcMm) 



(5) 



with k — 1, N/2 and N. All the averages are performed 
over 20 independent runs. 

We consider, first, the EPS for three values of the tilt 
F, first column of Fig. [31 For all values of F, the central 
particle moves with velocity very close to the center of 
mass velocity vcm as shown by the fact that stays 
almost pinned around the value 1. The velocities of the 
extremal particles, instead, are very different from vcai 



only during an initial transient when the system remains 
compact. Indeed, it is clear that, in a system not yet 
too diluted, the motion of the first particle is frequently 
hindered by the others, while the motion of the last is 
favored. Asymptotically, the differences in the motion 
of rods 1, N/2 and N become less evident because the 
interactions become less effective as the packet spreads 
over and over. Notice, also, the asymmetry between the 
first and the last particle. 

The major change observed when the inelasticity is 
turned on (a < 1) is that the velocity of the last particle 
displays a pronounced deviation from vcai- These dif- 
ferent behaviors are also evident by an inspection of the 
shapes of the instantaneous coarse grained distributions 
of particle positions (Fig. [4]), 



N{x,t) 



x+w/2 



-w/2 



dyp{y, t) 



(6) 



computed, in the simulations, by binning the number of 
particles to widths of the size of a single potential well. 

In the EPS, the asymmetry of the packet is determined 
only by the excluded volume effects and, when the diam- 
eter of the rods increases from a to lOcr, both asymmetry 
and drift velocity become more pronounced. 

The dissipation, favoring clustering, makes the packets 
of the IPS more compact and, in addition, changes the 
direction of the asymmetry, as seen for a = 0.9 and a — 
0.8, in Fig. [D This effect can be quantified by means 
of the skewness of the particle distribution around their 
center of mass. 



Skew(a;i, xjv) 



1 ^ 



rmsd 



XCM ^ ' 



(7) 



where rmsd(a;i, a;Ar) is the standard deviation and 
xcM the center of mass position. Positive values of Skew 
entail distributions with an asymmetric shape extending 
out toward more positive tails. Negative values are as- 
sociated with distributions extending out toward more 
negative tails^" . In the inset of Fig. [4l we plot the pa- 
rameter Skew as a function of time for three systems 
with different inelasticities. 

Let us consider, now, the quantity 



D 



lim -R{t) 

t — ^00 zt 



(8) 



In the IBP, D is constant and corresponds to the diffu- 
sion coefficient. We have found numerically that such a 
behavior persists both in the EPS and IPS for all val- 
ues of a we explored, as clearly indicated by the linear 
growth of R{t) in the inset of Fig. [5l 

The interactions change quantitatively the dependence 
of D on F, in fact, as F varies, the diffusion coefficient 
displays a maximum for F F3, (see Fig. [5]), with a 
behavior similar to that found by Reimann et al. in the 
IBP—. As a decreases, however, we observe some dif- 
ferences. In particular, the larger diffusivity of the EPS 
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FIG. 3: Evolution of the normalized displacements from the initial position (see Eq. (0) of the trajectories referring to tagged 
particles k = 1 (full line), k = 128 (dot-dashed) and k = 256 (dots). The normalization is performed with respect to the 
displacement of the center of mass from its initial coordinate. The panels al,a2 and a3 refer to rescaled forces FjF-i ~ 0.71, 
bl,b2 and bS to F jF-i ~ 0.88 and finally cl,c2 and c3 to FjF-i, ~ 1.06. The system contains = 256 particles, the coefficients 
of the restitution are a = 1.0 (in plots al, bl, cl), a = 0.9 (in plots a2, b2, c2) and a = 0.8 (in plots aS, b3, c3). Data are 
averaged over 10 independent realizations and the remaining parameters are the same as in Fig. [1] 



with respect to the IBP is basically explained by the ex- 
cluded volume repulsion. At first, the diffusivity lowers 
as a decreases, but a further decrease of a determines a 
new enhancement of D (Fig. [5]) . Such a feature is a con- 
sequence both of the larger broadening and of the larger 
asymmetry of the packet occurring as the particle size 
grows (see Fig. |4]); as a matter of fact, there is a huge 
increase of the diffusion since, for F/F3 ~ 0.88, we found 
-D/Dq ~ 29.6 (these data are not shown). 

IV. UNDERDAMPED REGIME 

When the damping force, in Eq. ([2]), is smaller or com- 
parable with respect to the inertial term mx the situation 
changes drastically. 

At low friction, the particles can travel across several 
wells before being trapped in a minimum. Again, as a 
guide, we can use the results obtained by Risken in his 
thorough study on the IBP in washboard potentials, pro- 
viding a full theoretical treatment of the mobility as a 
function of 7 and temperature. Risken's theory shows 
the existence of two dynamical states: locked and run- 
ning. At T = 0, a characterization of the dynamical be- 
havior of the non interacting system is straightforward. 
Indeed, for F > ^3 = only running states exist. When 
Fi < F < F3, with Fi — AjTr^/niVo particles may be ei- 
ther locked or running depending on the initial velocity 
and position. This fact determines a hysteresis loop in 



a fi vs F diagram. When the temperature is finite, the 
particles can switch from one state to the other under 
the influence of the thermal bath. As a result, the hys- 
teresis is suppressed and the locked- running transition 
occurs smoothly as a function of the tilt F, being the 
/i(T) curve a sigmoid in F. Only as T — > 0^ the sigmoid 
becomes a step function whose discontinuity is located at 
F2 ~ S.SGj^/mVo. 

MD simulations of the underdamped regime at differ- 
ent values of F allowed us to determine the curve of fi 
as a function of the tilt, shown in Fig. [6] for the EPS 
(a = 1) and the IPS with a — 0.8 at temperatures T — 8 
and T = 4. 

The mobility, fj,, appears to be roughly similar to that 
of the overdamped regime, but the critical F separating 
locked and unlocked situations lies at values lower than 
F3 and depends both on 7 and T. In the EPS, the mo- 
bility is rather close to the IBP value, while for the IPS 
the mobility is reduced, as shown explicitly in the case 
a = 0.8 in Fig. El 

While, the mobility of particles interacting inelastically 
displays, in the underdamped regime, no peculiar behav- 
ior with respect to the IBP and EPS, their collective dif- 
fusion presents some anomalies. At low temperatures, 
the dynamics strongly depends on the initial conditions, 
the evolution of the packet is very sensitive to the tilt and 
exhibits a rather inhomogeneous and irregular structure. 
This can be readily visualized by looking at the time 
behavior of the coarse grained particle density N{x,t) 
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FIG. 4: Averaged number of particles, for different coefficients 
of restitution, as a function the dimensionless quantity xj-w 
with F jFz 0.88 and — 2 ■ 10*. The circles correspond to 
the elastic case (a = 1.0) with a = 1, while the squares and 
triangles to a = 0.9 and a — 0.8, respectively. The diamonds 
correspond, instead, to the elastic particles of size lOcr. The 
distributions extend only over few hundred wells, while the 
system size is much larger L = lO^w. Inset: The trend of 
the skewness as a function of 'yt in the same cases. The other 
parameters are the same as in Fig. [T] 




FIG. 5: The rescaled coefficient of diffusion D/Do, {Do = 
T/mj) as a function of F/F3. Each point is the result of an 
average over 10 independent runs. The circles correspond to 
the elastic case (q = 1.0), while the squares and triangles to 
a = 0.9 and a = 0.8, respectively. The full line represents 
the theoretical result obtained by Reimann et al.— . Inset: 
standard diffusion of the system rescaled spreading, R{t)/a'^, 
[see Eq. ((T}] with respect to the dimensionless time 'yt in the 
same cases for F/F3 ~ 0.88. The other parameters are the 
same as in Fig. [1] 



[Eq. ©J at temperature T = 4.0 shown in figures [7]^a) 
and [Zl^b) for F = 1.4 and F = 1.6, respectively. It indi- 
cates that, for F = 1.4, particles remain partially trapped 
in the well where they have been initially deposited and 
only a fraction of them escape acquiring a drift. For 
comparison we also show the corresponding behavior of 
N{x,t) for the EPS, where the diffusion is standard. In- 
terestingly, the inelastic distribution N{x,t) corresponds 
to a lower mobility, but to a larger spread of the parti- 
cles. Moreover, one can see that new clusters, indicated 
by the spikes in Fig.[7{a), spontaneously form and persist 
for long periods before being dissolved. Such clustering 
phenomena are favored by the moderate value of the tilt. 
The situation, instead, looks different at F = 1.6 [see 
Fig-EIb)], where the initial cluster "evaporates" earlier, 
and the formation of new clusters is prevented by the ef- 
fects of the drift F. One observes a clear difference with 
respect to the behavior of the overdamped IPS illustrated 
in Fig. m showing a more compact structure of N{x, t) 
which does not lose "debris" . 

Such an early stage is sufficient to determine a late col- 
lective transport characterized by an anomalous spread- 
ing. The closely packed initial configuration has deep 
repercussions on the late spreading, R(t). This quan- 
tity, despite a very long simulation, does not appear to 
reach a linear or any other simple functional asymptotic 
dependence on time (Fig. [8]). The situation described 
above corresponds to an early very steep growth of R{i), 
growth that becomes slower after a characteristic time 
T whose duration decreases with the tilt amplitude. As 
shown in Fig. [51 the vertical position of the knee strongly 
depends on the number of particles. By no means the 



R{t) shows the expected linear behavior of standard dif- 
fusion. Such a feature is present only in the inelastic 
systems. However, when the temperature is raised, even 
the IPS recovers a linear behavior as seen in Fig. [Sl^b). 
The estimated diffusion coefficient is shown in Fig. [HI 

The peak of D relative to the IPS is higher than the 
corresponding peak of the EPS. The larger value of 
the IPS peak can be associated with the shape of the 
phase-space distribution, P(xr,Vr), which is defined as 
the probability of finding a particle at distance a;,- from 
the center of mass and with velocity Vr with respect to 
center of mass velocity vcm ■ 

The enhancement of the diffusion in the IBP is de- 
termined by the locked-running bistability, i.e. by the 
existence of two peaks at velocities v = and v — F/j 
in the velocity distribution. At the transition both pop- 
ulations are present and the center of mass velocity does 
not represent the most probable velocity in the system. 
Hence, it is clear that in this situation the variance of 
the velocity distribution can be very large. On the other 
hand, when one of the two populations becomes domi- 
nant the variance tends to decrease. The spatial part of 
the distribution P{xr,Vr) does not play any role in the 
IBP as the x-dependence of the distribution is symmet- 
ric with respect to the center of mass coordinate. The 
scenario changes in the interacting cases, since in one di- 
mension, the excluded volume plays a fundamental role 
in the dynamics. As we have seen in the inelastic case, 
at low temperatures the interaction was such as to lead 
to a fragmentation of the system and to a non linear be- 
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FIG. 6; Dimensionless mobility 7/^ as a function of the ex- 
ternal rescaled force F/F2 for different temperatures and co- 
efficients of restitution. The open and closed symbols refer 
to T = 4 and T — 8, respectively. The squares indicate the 
elastic case and the triangles the case with a = 0.8. The 
other parameters of the system are 7 — 0.2t^^, N — 256, 
Vo = 9.0, a = 1.0, w/a = 400 and t^ax = 10000t„. The full 
hne indicates the corresponding curve for the non interact- 
ing particles, obtained via the continued fraction method by 
Risken. 
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FIG. 7; Snapshots, taken a different times, of the averaged 
distribution N{x, t) plotted as a function of the dimensionless 
quantity x/w, for F — 1.4 (panel a) and F = 1.6 (panel b), for 
the IPS with a = 0.8 and EPS (broken line). The remaining 
parameters are the same as in Fig. [6] We notice that the 
system size is L = 10*?;; much larger than the well width, 
thus the peaks are not commensurate with the system size. 
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FIG. 8: Time behavior of R{t)/a^, 0, for the IPS with a = 
0.8. Panel (a) reports three runs with F = 1.4, T — 4 and 
— 256 (circles) two runs with A'^ = 128 (squares) and two 
runs with A'' = 64 (triangles). The full line is the average of 
the same quantity over 10 independent runs involving N = 
256 particles. Panel (b) shows R{t)/a^ vs -yt at temperature 
T — S for different forces: F = 0.8 (circles), F = 1.2 (squares) 
and F — 2.0 (triangles). The curves represent an average over 
five runs. The remaining parameters are the same as in Fig[6] 
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FIG. 9: Rescaled diffusion coefficient D/Dq {Do = T/m-f), 
for different coefficients of restitution, as a function of F/F'z. 
The D values represent averages over ten trajectories. The 
squares correspond to the elastic case (a = 1.0), while the 
triangles to q = 0.8. The temperature is T = 8.0, the re- 
maining parameters are the same as in Fig. [S] The dashed 
line indicates the result relative to the IBP at a lower temper- 
ature T — 4.0, showing that the effect of inelastic collisions 
cannot be accounted for by a mere reduction of the effective 
temperature of the system. The inset, finally, displays the 
ratio between the kinetic and the heat bath temperatures. 



havior of the diffusion. As T increases the fragmentation 
decreases and R{t) display a linear behavior with respect 
to t. In order to explain the different values of D near 
the peak between the IPS and the EPS, we have collected 
from the simulations the double histograms P{Xr,Vr). 
The comparison is reported in Fig. [TO] for drift values 
F near the peak. Analyzing the different behaviors of 



P{xr,Vr) in the various cases may help clarify the role 
of the interactions. In the EPS for loads F ~ 1.2 and 
1.4, P{xr,Vr) displays the same velocity bistability as 
the IBP and no spatial asymmetry and indeed the diffu- 
sion in these two systems results very similar as shown 
in Fig. [51 In the IPS, instead, for the same loads we 
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observe a strong spatial asymmetry which enhances the 
diffusion. In the IPS, the interactions besides determin- 
ing the peak of the diffusion also determine its width. 
From our simulations, as shown in Fig. [HI we see that the 
peak occurs at lower values of F with respect to the IBP 
and EPS and the width is much broader. This means 
that it is not possible to reproduce the behavior of the 
IPS by an appropriate choice of an effective temperature 
accounting for the inelasticity. A smaller temperature, in 
fact, would give a higher peak, but would not change the 
value of F at which it occurs. Moreover, the width of the 
peak decreases with T and fails to reproduce the observed 
broadening in the IPS, which is due to the inter-particle 
interactions. In the inset of Fig. [9l we display the ratio 
between the kinetic temperature 
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and the heat bath temperature T. Notice that in both 
cases the peak occurs in correspondence of the peak of 
the curve D/Dq of the elastic system and the ratio ap- 
proaches, as expected, the value 1 at low and high bias 
F. 
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FIG. 10: Joint distribution P{xr,Vr) of position and velocity 
in the reference frame of system center of mass, plots refer to 
force F = 1.2 (left) and F = 1.4 (right) in the EPS case (top) 
and IPS case (bottom) with a — 0.8. Independent variables 
have been made dimensionless through the rescaling x/xcm 
and v/vth, with Vth ~ \/T jm, the thermal velocity. 



mobility and the spreading of the system with respect 
to its center of mass, indicators that are customarily 
used to characterize collective transport is such systems. 
We compared the behavior of three models involving re- 
spectively: independent Brownian particles (IBP), elas- 
tic particles (EPS) and inelastic particles (IPS). We have 
found that, in the high friction regime, the transport be- 
havior of the system with interacting particles is qualita- 
tively similar to that of the non interacting system (IBP). 
Basically, the interactions do not strongly affect the mo- 
bility which remains close to the IBP value, while they 
modify the diffusivity that appears to be larger. This 
conclusion is intuitive, if we consider that p. is related 
the center mass velocity of the system, a quantity rather 
insensitive to the presence of interparticlc interactions. 
To achieve a a finer information about the structure of 
the collective motion, we monitored the evolution of the 
particle spatial distribution. Such a distribution strongly 
depends on the choice of the geometrical parameters and 
of inelasticity. This analysis shows that an initial local- 
ized packet of particle spreads differently in the inelastic 
system from the elastic one. 

The scenario in the underdamped regime is more com- 
plex. The EPS systems, in fact, are similar to IBP either 
analyzing the mobility, or the diffusion of the particles. 
The inelastic interactions practically do not affect the 
mobility, whereas determine a marked change in the be- 
havior of the diffusion at low temperatures and forces. 
The spreading, R{£) (Eq. ([T])), presents a clear two stage 
behavior: it increases nonlinearly in a transient regime 
which is then followed by a nonlinear growth, as shown 
in Fig. [51 At higher temperatures, however, the diffusion 
of the inelastic particles exhibits a single regime, with a 
linear growth typical of a standard diffusion. Our study 
indicates, moreover, that a group of particles interacting 
via inelastic hard-core collisions spread over much more 
than elastic particles. 

The mobility /i seems to be a global indicator that 
poorly encodes the information about the details of the 
interactions. On the contrary, the diffusion coefhcient D 
is a much more representative observable being more sen- 
sitive to the presence of elastic or inelastic interactions. 

We believe that in the overdamped regime, it should be 
possible to apply theoretical methods already employed 
in the study of dense molecular fluid in narrow chan- 
nels^, whereas the treatment of the underdamped case 
remains more problematic. The aim of the present nu- 
merical work is, perhaps, to stimulate further investiga- 
tion and new theoretical proposals. 



CONCLUSIONS 



In summary, in this work we have reported a numerical 
study on the dynamics and the transport of an array of 
identical particles in a inclined washboard potential. We 
investigated, at several tilting F, the role of hard-core 
repulsive interactions among the particles on both the 
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